Generation of pure bulk valley current in graphene 
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The generation of valley current is a fundamental goal in graphene valleytronics but no practical 
ways of its realization are known yet. We propose a workable scheme for the generation of bulk valley 
current in a graphene mechanical resonator through adiabatic cyclic deformations of the strains and 
chemical potential in the suspended region. The accompanied strain gauge fields can break the 
spatial mirror symmetry of the problem within each of the two inequivalent valleys, leading to 
a finite valley current due to quantum pumping. An all-electrical measurement configuration is 
designed to detect the novel state with pure bulk valley currents. 

PACS numbers: 72.80.Vp,85.85.+j,73.63.-b 



Apart from pseudospin (chirality) , charge carriers in 
graphene are also characterized by the valley index 
(sometimes called isospin) originated from the existence 
of two conical (Dirac) points per Brillouin zone[TJ. A val- 
ley polarized state requires the absence of time reversal 
symmetry, as the two valleys are related by this symme- 
try. 

Motivated by the growing field of spintronics[2], it 
was proposed that the manipulation with the valley 
index may open a new way to transmit information 
through graphene, and different manipulation schemes 
were proposed [3T4TU]. After initial enthusiastic attitude, 
interest in "valleytronics" declined somehow, as it was 
soon realized that a valley polarized current will be de- 
graded by intervalley scattering induced by atomic scale 
disorder lj, making it difficult the maintenance of val- 
ley polarized states. In addition, a number of proposals 
were based on the spatial separation of valley currents at 
zigzag edges [3], which requires well-defined edge orienta- 
tion and at the same time free of short-range scattering. 

We present a new scheme to induce valley polarized 
currents in graphene which avoids some of the pitfalls of 
previous proposals. The breaking of time reversal sym- 
metry is achieved by means of time dependent fields, 
instead of a magnetic field. The induction of valley 
polarization by a.c. fields has been proven [TTj in M0S2 
where optical radiation was used in order to excite val- 
ley polarized charge carriers. As in this experiment, 
the scheme discussed below generates valley currents 
throughout the entire system. However, the a.c. driv- 
ing force in our proposal is due to mechanical vibrations 
of a nanoelectrical- mechanical system (NEMSfT^HH]), as 
illustrated in Fig. la. 

We employ a continuum-medium description of 
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FIG. 1: Illustration and definitions of the (a) graphene based 
NEM device and (b) crystallographic coordinate systems used 
in the paper. 



graphene with the Dirac Hamiltonian, 

•H ± (fc,^') = hv f [±(k x TA x )a x + (k v TAy)a y }+fit 2 , (1) 

where +(— ) denotes K(K') valley index, Vf is the Fermi 
velocity and \x is the chemical potential. A{f) is an ef- 
fective gauge field describing the modifications to the 
hopping amplitudes induced by the strain fields u^{f) 
[TBI H7] , and has opposite signs at the two valleys as re- 
quired by time-reversal symmetry. We described the de- 
formation of the suspended region, i.e. —U<y<^, 
with a simple uniaxial strain given by u yy — u and 
u xx = u xy = 0. It will be shown below that for arbi- 
trary crystallographic orientation 6 (see Fig. lb), A(r) in 
the suspended region is given by the expression 



XCf) = ^(-cos30,sin30), 
a 



(2) 



where ft ~ 2 is the electron Gruneisen parameter and 
parameter related to graphene's elastic prop- 
erty as described in [16] . a ~ lAA is the interatomic 
distance. Adiabatic cyclic variations of the internal pa- 
rameters, such as deformations in the strains (it) and 



2 



chemical potential (fx) in the suspended region, over a 
work cycle can constitute a scheme for adiabatic quan- 
tum pumping [I5|IT5] (for the general theory of quantum 
pumping, see Ref. [201 [21] )• For the charge pumping 
one needs to break the inversion symmetry, e.g., mak- 
ing the right and left leads different (e.g., by different 
doping) [21 . Here we will demonstrate that for the case 
of symmetric leads the valley current through the system 
is, in general, nonzero. In this case, the pumping cur- 
rent through each channel will be shown to follow the 
general relation (with periodic boundary condition along 
the transversal x direction), 

Il^ K (k x ,6) = -Il^ K '(-k x ,9), (3) 

where the subscript L/R refer to left/right leads(/ £ /^ is 
defined postive when the current flows out of the device). 
This significant relation results from an intrinsic sym- 
mety of graphene NEMs, as we will show later. Eq.([3| 
represents the central result of our paper. Below, we 
present the derivations leading to Eq.([3]) and discuss the 
physical consequences that follows, such as the genera- 
tion of pure bulk valley current and its possible experi- 
mental detection. 

We start by considering the case 9 = 0, where the 
trench is directed along the zigzag direction (see the 
(£o> Vo) coordinate system shown in FigflJ)). In this case, 
the Hamiltonian has the form of Eq.Q and the pseudo- 



where we have defined the pseudo-magnetic field A in the 
rotated frame as |22| : 

{X) =n{W) {^^~ul VV) )- (8) 

Eqs.Q and ^ constitute the continuum-medium de- 
scription of stained-graphene in an arbitrarily rotated 
frame. Here follows two comments about this general 
form. First, if we perform the gauge transformation of 
the wave function on B sublattice ipg K — > e ±rf , 
the Dirac Hamiltonian can be made explicitly invariant 
(i.e. Eq.Q) under rotation. Thus, we can simply drop 
the factor e ±l ° in Eq.(|7j) in subsequent discussion. Sec- 
ond, from the definition for pseudo-magnetic field, it is 
obvious that the form is of 27r/3 periodic in 9, which 
reflects the trigonal symmetry of the underlying honey- 
comb lattice. 

Next, we describe the quantum pumping problem 
based on graphene NEMsfTT)! |2"3"] . As discussed above, 



magnetic vector potential reads: 

A - - \ 

J\ x — \U XX W>yy ) i 

2(3 K 

•Ay — U xy . (4) 

For arbitrary orientation i.e. 9^0, Eq.Q and Eq.Q 
have to be recasted in the new coordinate frame (x,y). 
The two coordinate system are related by 

( x )=n( x °) , n=( cos ° sin jV (5) 
V y J \ya J V - sm # cos # / 

The wave vector k transforms in the same way as r such 
that k ■ f is a rotational invariant quantity. The strain 
field is defined as, Uij = + ff 1 ), which is a sym- 

metric tensor of rank two. Hereafter we use the sub- 
script/superscript "0" to denote physical quantities in 
the original frame (xo,yo). Thus, we have ko = R~ 1 k 
and 

u xx - u yy = cos 29(uxx - u yy ) - 2 sin 29u xy , 

v° xy = - sin 29(u xx - u yy ) + cos 29u xy . (6) 

By using the new coordinates in the Dirac Hamiltonian, 
we can transform it to the rotated frame: 



(7) 

I 

the derived Hamiltonian given by Eq.([7]) is physically 
equivalent to Eq.|lj). Suppose that between — ^ < 
y < U, the system has uniaxial strain u yy = u and 
Uxx — u xy = 0, hence, Eq.Q is reduced to 

(A x ,A y m = I caafywiM), \y\ < J (9) 

with the expression given in Eq. ^ for the suspended re- 
gion. In the following, assuming a particular geometry 
(the width W^>L), the x direction is treated as trans- 
lationally invariant. From Eq.([7| and Eq.Q, one can 
easily see that 

H + (k x ,k y ,A,9) = H-(-k x ,k y ,A,-9). (10) 

We can call such combined symmetry as the crystalline- 
angle- combined mirror symmetry in the continuum- 
medium model. It turns out that such combined sym- 
metry has a significant consequence on the relation of 



l-L ± (k.A) = hv f 



,±i6 



[(±k x + ik y ) T (±A X + iA y )\) 



'[{±k x - iky) =F (±Ax - iA y )] 
fi/Hvf 
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pumping currents in each valley, i.e., Eq.([3]), as will be 
elaborated further. 

Due to the mentioned symmetry, we may focus only on 
the K Dirac cone, whose Hamiltonian has the following 
form for the different regions (i.e., i — L,R,G denotes 
y < -k , y > § , \y\ < § , respectively): 



'(k,A) = hv f ( ,~ ei/hv J . (fc * /b lky)t ) , (11) 
T \(k x + ik y )i a/hvf J ,K ' 



where we have defined (k x , k y )i = (k x — A x , k y — A y )i- 
In this paper, we consider mainly the symmetric case 

The eigenenergies and eigenstates of the Hamiltonian 



in Eq.( 11 ) read 



Ei(k) 



< ; + shVf\l kl + k%, 



(12) 



where s — ±1 refers to electron/hole band, respectively. 
Due to translational invariance in x direction, k x is the 
same in all three regions. We consider now the equilib- 
rium situation where all three regions can be described by 
a common chemical potential /i. Obviously, k Xi is always 

real since k 



is real. Then, k 



± 



k 2 ., where 



kf t = j-^- and the ± sign is selected to give the correct 
sign in the group velocity, depending whether it is an 
incident, transmitted, or reflected waves. Note that k y 
can be purely imaginary representing evanescent waves 
in the central region. 

It is straightforward to calculate the scattering matrix 
for our device. Without loss of generality, we can focus 
on the case with the electron doping (fi— cl > 0) in the 
leads. The scattering coefficients are calculated to be 



£ -iky L L 



C 2 (d) + C 3 (d) 
C\(d) ' 
-ik VL L "4 sin ip L sintpge 

Ci(d) 



A-uLd 



(13) 



where and t& are the reflection and transmition coef- 
ficients with d = 1(— 1) corresponding to the cases with 
incident waves from y = — oo (left lead) and y = oo (right 
lead), respectively; ipi, and tpc arc defined through 



e ±iVL = 



k x ± ik VL d ±i VG _ k XG ± i 
toft, kfc 
defined as 



krr^ ± ik ya ^. 



The Ci(d) in Eq.([l3) are 

C\{d) — 4i sm(k VG Ld) (1 — cos ipL cos pc) 



(15) 



\--y G 1 \- 7- jj 

4cos(k VG Ld) sin cpi, sin ipo, 
-2i(l + e 2lipL )sm(k yG Ld), 

Ld-ip G )] 



C 3 (d) = -2z(l + e : 

C 3 (d) = 2ie iVL [sm(k VG Ld + if G )+sm(k yG 



Symmetries related with td and r^. Next, we discuss 
symmetry properties of the scattering amplitudes. First, 
we note that <PG(d) = <£>g(— G0, <Ph{d) — —(pL{~d) and 
they are independent of the sign of 9. Then, we can 
obtain the relations satisfied by CVs: C\(d) = —C\{—d), 
C*2,3(— d)* = C2fi{d). Based on these relations and the 
odd parity of A y (9), we arrive at: 



t d (k x ,0) 
rd{k Xl 9) 



t-d(k x , —9), 
rd{k x , -9). 



(16) 



Symmetry of pumped valley- dependent current. Ac- 
cording to the adiabatic pumping theory [20] and the sym- 
metry relations satisfied by and td, we obtain the fol- 
lowing relation for the pumping current I^ mp ' K {k x ,9) 
for K valley: 



jpump ,K 



ieui 
4^ 



jpump,K 



2 ^ dn{6) 

ds[ ^s- ri{d) 



(k x ,-9) 



dr-!(0) 
ds 



r-m% (17) 



where we've used s as time symbol to avoid confusion 
with transmission coefficient. Now, using the symme- 
try relations satisfied by CVs and the fact that the com- 
mon complex factors for Ci and C3, i.e., ie lipL , is in- 



dependent of time(see Eq.(19l for typical time depen- 



dence of pumping parameters for graphene NEM), we 



can simply prove that the integrand in Eq.(17l is zero 



By further taking into account the current conserva- 



tion condition( I 



pump,K 



(k x ,-6) 



0) and the symmetry relation guaranteed by Eq.(10) 
(i.e., jP ump ' K (—^0) — jpump.K —8)), we arrive at 
Eq.(|3]), the relation for pumped valley-dependent cur- 
rent, which is the central result of this paper. 

Next, we discuss some direct consequences of Eq.Q. 
By integrating out the transversal wave vector k Xl we can 
get the following relation, 



jpump,K' ,n\ 



(18) 



Eq.( 18 ) means that the total pumping current at a given 
lead is opposite for different valleys. Thus, the total 
charge current is exactly zero. This situation is analogous 
to the pure spin current generation in spintronics |241 [2"5] . 
thus we call this effect pure valley current generation. 
In summary, we have shown that the application of an 
alternating back gate voltage to graphene NEMs can in- 
duce a pure valley current via adiabatic pumping. In 
adiabatic pumping theory, there are two necessary con- 
ditions for finite charge pumping effect, i.e., time rever- 
sal symmetry breaking and mirror symmetry breaking 
of the the left/right leads [21]. The above derivations 
show that the valley pumping effect can be realized in a 
seemingly symmetric two-terminal geometry with identi- 
cal leads. However, the mirror symmetry of the system 
(i.e. y^—y) has actually been broken for each valley 
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FIG. 3: (Color online) Pictorial sketch of the pumping gen- 
eration of pure valley current and an all-electrical detection 
scheme. The Hall voltage difference V12 = Vi — V2 and 
V34 = V3 — Vi across the NEMs is predicted to bear oppo- 
site sign due to the flow of pure valley current. 



FIG. 2: (Color online) The pumped charge per cycle of the 
K valley, denoted as ■ The valley current can be obtained 
by eujQ* /2it. The width is fixed to be 5000a (w 700nm), 
calculated for two different length as indicated. We fixed the 
phase difference of the driving parameters <j) = tt/2 and take 
I j4q.ii I /A: / ~ 0.002, which amounts to a strain u ~ 1.75 x 10~ 4 . 



nent. For comparison, we have jP™ m P' K / K ( 
when A y = 



when Such symmetry breaking is embodied 

in the continuum theory through a nonzero A y compo- 

= 0) = 

which can be infered from the relation 

I p L ump - K/K '(k x ,9) = -P™ mp < K/K '(k x ,-6) implicit in the 
discussion below Eq.(17). As is worthy to recapitu- 
late, the symmetry property intrinsic to the graphene 
NEM(i.e., the crystalline-angle-combined mirror symme- 
try) can impose strong restriction of pumping currents 
on different valleys. For more quantitative understand- 
ing, we present some numerical results of jP^ m P' K ($) ( m 
terms of the pumped charge per cycle ) for the K valley 
in Figj2] using some typical experimental parameters, as 
discussed below. 

As stated above, the strain (u) and Dirac potential 
(ec) in the suspended region are modulated by the ac 
back gate voltage. Near resonance, they differs by a phase 
difference <fi with a typical time dependence given by(with 
conventional time symbol t) [H] : 



£g = £Go[l + 6ecos(u;i)]2, 
A = A [l + 6Acos(uit + (f))] 2 . 



(19) 



Assuming typical numbers for the static part eco — 
(-L = 0.3eV, Ao — 0.02kf(cos39, — sin 30), and the os- 
cillating amplitude Se = SA = 0.2 , we calculated the 
pumped charge per cycle (9) for K valley, as shown in 
Fig. 2. By definition, the pumping current jP^ m P' K ($) = 
euQ^ (8) /2tt. Our calculation indicates a nearly linear 
scaling(unshown) of the pumping effect on length L of 
the NEMs, which is similar to results infTU]. As explic- 
itly shown, the maximum pumping effect is reached for 
the crystallographic angles corresponding to armchair- 



type x axis {9 = tt/2 + rnr/3, n € Z) while it is zero at 
zigzag-type x axis {9 — nir/3,n G Z). The periodicity 
27r/3 with 9 is easily seen. The valley current can be 
defined as lP um P-^(9) = IP um P> K (9) - IP um P- K ' (9), which 
is simply twice the value of jP um P< K (9). For typical reso- 
nance frequency of w w 10MHz[T3] to 0.16GHz [TS], we 
arrive at numerical estimates IP um P' v « 0.1 — 10pA//im, 
a quantity measurable in experiment. 

The possibility of pure bulk valley current in this sim- 
ple pumping scenario looks very promising. The prob- 
lem is how to probe the valley current. Here we propose 
an all-electrical measurement as shown in Fig(3j From 
we can infer that the charge current pertaining 
to carriers from the valley K not only has opposite lon- 
gitudinal component with respect to the charge carriers 
from the valley K' , but also their transversal velocities 
are opposite. As pictorially shown in Fig(3j we expect 
charges accumulating on opposite edges on the two sides 
of the NEMs. Based on this observation, we predict that 
the resultant Hall voltage on the left lead has opposite 
sign with that on the right lead, i.e., Sign(Vi2/V34) = — 1, 
which is the characteristic feature of the bulk valley cur- 
rent flow. 

The above discussion can be made more clear if we con- 
sider the valley-dependent pumping Hall current, which 
can be calculated as l"^ l / K {9) = £ fex 7™'* 0)j£. 
From Eq.(|3]) and the particle conservation law, we can 
obtain: 



T Hall,K { r.\ 



T Hall.K / f\\ 
l LjR \°> 



T Hall.K / r\\ 
- J R/L (°>- 



(20) 



Such Hall current on different leads can result in the op- 
posite Hall voltage due to the charge accumulation near 
edges. Thus the existence of valley polarized currents 
possibly can be detected in the leads by means of nonlo- 
cal multiterminal measurements [26] . 

Our discussion above is restricted to the case with sym- 
metrical leads. It is straightforward to extend our study 
to the more general case with differently doped leads. 
In the general situation, the current is not purely valley 



5 



current, i.e., I P ™ P ' K {6) ^ (0), thus the charge 

pumping current is finite. An example is again given by 
the case 6 = 0, where I pu "^ K (0) = I%% P>K ' (0) M and 
the valley current is zero. 

To conclude, we have shown that through pumping 
induced by mechanical vibrations bulk valley polarized 
currents can be generated in graphene leads connecting 
the graphene resonator with trench directed at a gen- 
eral crystallographic angle. We have demonstrated that 
the generated current is purely valley current (with zero 
net charge pumping current) in the setup with the same 
doping rate in the graphene leads. Together with an all- 
electrical measurement scheme, our proposal opens a new 
direction of exploiting the valley degree of freedom, thus 
pushing forward graphene-based valleytronics a step for- 
ward toward real applications. 
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